Dynamics of a domain wall and spin-wave excitations driven by a mesoscopic current 
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The dynamics of a domain wall driven by a spin-polarized current in a mesoscopic system is studied 
numerically. Spin-mixing in the states of the conduction electrons is fully taken into account. When 
the Fermi energy of the electrons is larger than the exchange energy (Ef > J s d), the spin precession 
induces spin wave excitations in the local spins which contribute towards the displacement of the 
domain wall. The resulting average velocity is found to be much smaller than the one obtained 
in the adiabatic limit. For Ef < J s d, the results are consistent with the adiabatic approximation 
except for the region below the critical current where a residual domain wall velocity is found. 

PACS numbers: 75.75.+a, 72.25.Pn 



X 



The interplay between the dynamics of conduction 
electrons and the local spins in magnetic multi-layers has 
been pioneered by Berger 0,0] and Slonczewski |j. They 
pointed out that spin flips of conduction electrons induce 
a spin torque acting on the interfaces between multilay- 
ers. This spin torque can excite spin waves propagating 
in the interface perpendicular to the current direction. 
In the case of a smooth domain wall, as a result of a cur- 
rent flow, the spin torque can lead to domain wall motion 
opposite to the current direction. The spin torque can 
also induce an out-of-plane component of the local spins 
if damping is present. In addition, the reflection of the 
conduction electrons can lead to a reaction force acting 
on a domain wall via a momentum transfer. This can 
also lead to a displacement 0. 

During the past decade, considerable progresses have 
been made in understanding the current-induced domain 
wall motion in magnetic nanowires 0, 0, IE S 13 • 
This was motivated by remarkable experimental suc- 
cesses measuring the motion of a domain wall under the 
influence of a current pulse. Velocities of about 3m/s 
have been found which exhibit large fluctuations around 
the average value. The velocity of a domain 

wall induced only by the spin torque has been studied 
0, 0- This is suitable in the adiabatic limit, which is 
assumed to be a good approximation for ferromagnetic 
nanowires. In this limit, the thickness of the domain wall 
is assumed to be much larger than any of the character- 
istic length scales of the conduction electrons. The mo- 
mentum transfer is then very small due to the absence of 
backscattering. A velocity of 250 m/s has been estimated 
in lowest order of the interaction between the electrons 
and the wall. 

In the adiabatic limit, a finite final velocity of the do- 
main wall can only be sustained by the presence of an 
additional external magnetic field |4| and/or by lowest- 
order non-adiabatic corrections towards the shape of the 
magnetization profile 0, ||| . The effects of these correc- 
tions reduce the domain wall velocity at least by one or- 
der of magnitude. The current-driven dynamics of a do- 
main wall including contributions from non-adiabaticity, 



has been investigated microscopically by treating the 
conduction electrons quantum mechanically and includ- 
ing both, spin and momentum transfers Spin wave 
generation has been neglected. A parameter-dependent 
critical current has been found below which the domain 
wall is pinned, and which vanishes for an abrupt wall. 

In the present paper, we report results of numerical 
investigations of the motion of a domain wall driven by 
an electric current in the complete range of parameters. 
The conduction electrons are treated quantum mechani- 
cally within a lattice model. Quantum interference, spin 
mixing, and spin wave generation are taken into account. 
The time evolution of the local magnetization is described 
using the Landau-Lifshitz-Gilbcrt equations (LLG) in the 
presence of the magnetization of the current carrying 
state. The electric (spin-dependent) current in the pres- 
ence of the domain wall is determined as a function of 
time by using the quantum mechanical transfer matrix 
approach. When the Fermi energy E-p of the conduction 
electrons (measured from the band edge) is larger than 
the exchange energy J s d between their spins and the local 
spins of the domain wall, the state of conduction electrons 
is spin mixed. Then interference effects lead to spin pre- 
cession in the system of the conduction electrons. This 
excites spin waves in the system of the local spins which 
propagate in the direction of the current. The spin waves 
distort the shape of the domain wall from the saturation 
configuration of the adiabatic approximation. As a conse- 
quence, the domain wall can move beyond the saturation 
displacement. However, the resulting average velocity is 
much smaller than the one obtained from the spin torque 
only. In addition to this, we obtain strong fluctuations 
of the velocity which, to the best of our knowledge, have 
not been addressed before. For E-p < J s< \, the minority 
spin state becomes a damped mode outside the domain 
wall region. The spin precession of the conduction elec- 
trons does not occur, and the spin mixed state merely 
mimics an external magnetic field. This induces a small 
residual domain wall velocity below the critical depinning 
current predicted by taking into account only the spin 
torque mechanism. Our results strongly indicate that 
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the influence of spin waves cannot be neglected when in- 
terpreting experimental data. In particular, we suggest 
that the small velocity of the domain wall observed in 
experiments for currents well below the critical current 
results from the non-adiabatic contribution of the con- 
duction electron spin distribution. This could be viewed 
as the microscopic mechanism behind the non-adiabatic 
effective field that has been suggested previously 0, . 

We consider a one dimensional ferromagnetic spin 
chain in the x-direction (spins S.; at the sites i) and a fully 
polarized electron propagating in the +x-direction cou- 
pled via s-d exchange. For the magnetization Mj oc Si 
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where a the Gilbert damping parameter and M s the sat- 
uration magnetization. The effective magnetic field 



H effi = ffex {Mi-! + M i+ l 
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acts on the local spin at site i, with fl cx and fl s d the ex- 
change constant between the nearest neighbor spins and 
the s-d exchange constant, respectively, and Hk and M y i 
the anisotropy constant and the demagnetization field, 
respectively. The expectation value (s) i is the s-d ex- 
change field due to the polarized conduction electrons. 

In the adiabatic approximation, the conduction elec- 
trons are assumed to propagate without reflection and 
their spin directions follow the local spins. The corre- 
sponding effective magnetic field is oc M x (V • j e )M 
whereje is the unit vector in the direction of the cur- 
rent |4j. This is suitable for describing a ferromagnetic 
metal where the Fermi wavelength is much smaller than 
the system size and the width of the domain wall. 

However, for the mesoscopic quantum system consid- 
ered here, one must take into account interference and 
spin mixing of the states of the conduction electrons. For 
magnetic semiconductors, for example, the Fermi energy 
is of the order of a few meV and the phase coherent 
length can exceed several fmx. In order to treat con- 
duction electrons in such mesoscopic systems we have to 
use numerical techniques. Without interaction, the most 
simple lattice Hamiltonian (lattice parameter a = 1) is 
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with J sc j > and cf(ci) denoting the creation (annihila- 
tion) operator of electron at the i-th site. The hopping 
energy in the first term is assumed to be the energy unit, 
and hopping is restricted to nearest neighbors. The <x 
are the Pauli spin matrices. By using the transfer matrix 
method |ljj, [lj] one can calculate not only the transmis- 
sion and the reflection amplitudes of the domain wall but 



also the spin resolved current carrying states ip C cs- The 
expectation value of the s-d exchange field is 

(s) = <t/> C cs| IV'ccs) • (4) 

We determine the time evolution of the magnetization 
according to the LLG by taking into account this field. 
Both, the spin reaction torque and the momentum trans- 
fer are included. However, for the parameters used here, 
it turns out, that the momentum transfer is very small 
because the reflection probability of the system is less 
than 10 -2 . Since the current carrying state is influenced 
by the variation of the local spins, we have calculated the 
current carrying state in each time step. The calculation 
proceeds as follows. (1) An ad hoc initial (t = 0) spin 
configuration of the local spins Si of the domain wall is 
assumed. (2) The s-d exchange field (s)^ is calculated for 
a given magnetization by the transfer matrix method. (3) 
The time evolution of the local spins is calculated taking 
into account (s) i . (4) After a period t = 6t one obtains a 
new configuration of local spins Si from LLG. (5) Return 
to step (2) and calculate (s) i by using the new configura- 
tion. Here, we do not discuss the origin of the initial do- 
main wall. In principle it can be determined by a classical 
Monte-Carlo simulation for the given geometry of the sys- 
tem. The current, and thus the time dependent electric 
conductance contains important experimentally accessi- 
ble information on the dynamics of the domain wall which 
will be reported elsewhere fl5j . The initial spin configu- 
ration is assumed to be a smooth function of the position, 
S = (,S[cosh(X DW - x)/VF]- 1 ,0,S , tanh(A DW - x)/W), 
where W and A"dw are the width and the center of the 
domain wall, respectively. We have assumed W = 5, 
Adw = 30 and 5 = 1. The resulting s-d exchange field 
is shown in Fig. ^ 

We consider two energy regions: the adiabatic, E-p < 
J s d, and the mesoscopic region, Ep > J sc i- In the former, 
the wave number of the minority electron spins is purely 
imaginary outside the domain wall and the spin is almost 
parallel to the local domain wall spin. Therefore, the out- 
of-planc component (y-component) of the s-d exchange 
field exists only in the vicinity of the domain wall where 
the conduction spin rotates in the (a;-z)-planc. In the 
mesoscopic case, the spin-mixed state is allowed in entire 
system. Quantum interference of the spin states causes 
spin precession of the conduction electrons (Fig.[TJ. This 
leads to generation of spin waves in the system of the 
local spins. 

The time evolution of the domain wall is obtained by 
using a 4-th order Runge-Kutta method. The structure 
of the magnetization is shown for t = 100 in Fig. ^ 
The parameters are flex = 3.0, fl s d = 0.1, flu = 0.1, 
M s = 1.0 and a — 0.03. Fixed boundary conditions are 
assumed at the edges of the system of the local spins. One 
notes the displacement of the domain wall and that the 
out-of-planc component of the domain wall is developed 
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FIG. 1: Top: The initial s-d exchange field of the conduc- 
tion electrons, (s), calculated by the transfer matrix method. 
Bottom: The local spin configuration for t = (dashed dot- 
ted lines) and after time evolution (t — 100). Parameters: 
-Hex = 0.3, s-d exchange H s d = 0.1, anisotropy field Hk = 0.1; 
adiabatic region Ef — 0.11, J s( j = 0.2 (dotted lines); meso- 
scopic region Ef = 0.21, J sc i = 0.2, a = 0.03 (full lines) 



in both regions. In the adiabatic case, saturation of the 
out-of-plane component is obtained, in agreement with 
previous works In the mesoscopic case, one observes 
spin waves that have been excited via the s-d exchange 
field. These cause a distortion of the domain wall because 
they propagate along the current direction. 

Figure |21 shows the center of the domain wall as a func- 
tion of the time. We define the center of the domain 
wall as the position where M z = 0. In the adiabatic 
region, the displacement becomes saturated after a rel- 
atively short period of time. The saturation of the do- 
main wall motion is due to dissipation via damping. In 
the mesoscopic case, the displacement shows fluctuating 
behavior that corresponds to regularly accelerating and 



slowing down of the domain wall during the time evo- 
lution. This is due to the excitation of the spin waves. 
They induce distortions of the wall from the saturation 
configuration that is expected in the adiabatic limit. As 
a consequence, after the onset of saturation due to the 
damping, the wall again starts accelerating due to the 
distortion process before it gets slowed down again as 
a result of the damping, and so on and so forth. The 
high-frequency oscillations in the motion reflect standing 
wave patterns of spin waves induced by the boundary 
conditions. 




FIG. 2: The position of the center of the domain wall (Xdw) 
as a function of time for the adiabatic (dotted line) and the 
mesoscopic case (full line). Parameters: sd-exchange H B a = 
0.04, anisotropy field H K = 0.1, a = 0.03. 

Now we focus on the velocity of the domain wall. Fig- 
ure 01 shows the time average as a function of the s-d 
coupling constant H s d . The latter is proportional to the 
current. In the adiabatic case, we have plotted the ini- 
tial velocity because the motion saturates during time. 
It has been predicted that the velocity is proportional 
to \J H^ d — Hf t , where H CI is the critical s-d exchange 
constant which is proportional to the depinning current 
0- Our results are consistent with this for H s d > H CI . 
However, for H s( \ < H cv we obtain a small but finite resid- 
ual velocity as shown in Fig. In this region, the mo- 
tion of the domain wall saturates also, but after a much 
longer time than above the critical current. The resid- 
ual velocity is due to the non-adiabaticity induced by the 
non-adiabatic spin mixing in the state of the conduction 
electrons @. The s-d exchange field and the local spin 
are separated by a small angle inside the domain wall 
region. Thus, the domain wall feels an effective magnetic 
field that corresponds to the out-of-plane component of 
the s-d exchange field. The effective magnetic field in- 
duces a displacement of the wall even below H cr where 
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the spin torque is canceled out by the anisotropy field. 

For the mesoscopic case, we have plotted the average 
velocity and the square root of its variance because the 
propagation of the wall includes fluctuations, as shown 
in Fig |21 For H s( i > i? cr , the velocity is small compared 
with the initial velocity obtained in the adiabatic case. 
This discrepancy can be explained by considering the in- 
terplay between the spin torque and the spin waves as 
mentioned before. In the adiabatic limit the spin torque 
is the only mechanism that contributes. 
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FIG. 3: The domain wall velocity (V) as a function of the 
s-d exchange constant (i/ s d) that is related to the amplitude 
of the current. The initial velocity is plotted for the adiabatic 
case (Ef < J B a, Ef = 0.11, J s( j = 0.2, bullets and dotted 
line) while the average velocity is shown for the mesoscopic 
case (-Ef > J s d, Ef = 0.21, J s( j = 0.2, circles). Error bars: 
square root of the variance of the velocity. Solid line: the 
analytical result oc * / H^ d — H% T , where H CT is the critical s-d 
exchange constant that is estimated as H CI ~ 0.014. 



In conclusion, we have numerically studied the dynam- 
ics of the domain wall driven by a mesoscopic current. We 
have treated the conduction electrons quantum mechan- 
ically. Spin-mixed states are taken into account. In the 
adiabatic case, Ep < J sc i, the contribution of the spin 
mixed state is equivalent to an additional external mag- 
netic field which induces a residual velocity of the domain 
wall. In the mesoscopic region, Ep > J s( j, the spin preces- 
sion of the conduction electrons induces spin wave exci- 
tations in the system of the local spins. We have pointed 
out that the presence of these spin waves enhances the 
displacement. The resulting average velocity of the do- 
main wall is small compared with the one obtained from 
the adiabatic spin torque mechanism only. In the limit 
of small driving current the velocities of the adiabatic 
and the mesoscopic cases coincide within the error bars. 
Our results suggest that non-adiabatic effects, that can 
always be expected to be present in experiment, give rise 



to distortions of the local spin configuration and lead to 
a finite velocity of the wall even for currents below the 
critical value. We suggest that this is the microscopic 
mechanism that accounts for the non-adiabatic contri- 
bution towards the shape of the domain wall discussed 
earlier 0, Q . This could explain very well the results of 
the recent experiments where rather small domain wall 
velocities have been found for currents that are an order 
of magnitude below the critical current 0, ■ 

An alternative explanation of the experimentally ob- 
served velocities could be the presence of spin waves 
which could be due to temperature effects even for an 
adiabatic system. This would lead to much smaller ve- 
locities than predicted adiabatically (cf. Fig.[3J). Further- 
more, from our results, the spin-wave induced velocity is 
predicted to show large fluctuations that are consistent 
with the experimental observations [lfjj . 

Finally, we want to point out that our calculations 
have been done for a general abstract model. By suit- 
ably adjusting the parameters, and by generalizing the 
model Hamiltonian, for instance by including more than 
one conduction channel, our method may be able to 
contribute to better understanding a variety of differ- 
ent experimental situations ranging from ferromagnetic 
nanowires to magnetic semiconductor nanostructures. In 
addition, since we treat the conduction electrons fully 
quantum mechanically, mesoscopic interference effects 
in the current as a function of time can in principle 
be expected — especially for magnetic semiconductor 
nanowires 0, 0] — that could be signaled in the dy- 
namics of the domain wall. Whether or not this will be 
the case is an open question and the subject of future 
research. 
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